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I.  INTRODUCTION 


A.  THE  PROBLEM 

The  problem  considered  here  is  to  defend  a  multiple-ship  Battle  Group  (  BG  ) 
against  anti-ship  cruise  missile  (  ASCM  )  attacks.  For  this  reason,  a  multiple-ship  mis¬ 
sile  attack  model  has  been  developed  to  study  the  effect  of  Electronic  Warfare  (  EW  ) 
expendable  decoys,  which  may  increase  the  survivability  of  a  single  ship  or  a  battle 
group.  The  effectiveness  of  a  single  decoy,  or  of  a  group  of  identical  or  different  decoys, 
employed  during  an  ASCM  attack  on  a  BG  may  be  degraded  as  a  result  of  missile  di¬ 
version  onto  neighboring  own-force  or  allied  ships  in  the  BG. 

The  model  that  we  present  in  Chapter  III  is  intended  to  investigate  how  expendable 
decoys  may  increase  survivability.  The  emphasis  is  on  simplicity  consistent  with  repres¬ 
enting  the  multi-ship  situation  rather  than  on  complete  realism.  Individual  systems  may 
have  peculiarities  that  may  make  some  of  the  assumptions  in  the  model  invalid. 

In  this  study  it  is  assumed  that  the  BG  is  being  attacked  by  one  or  more  guided 

missiles.  For  each  individual  missile,  an  attack  is  assumed  to  consist  of  two  phases: 

1.  Target  acquisition  by  the  missile  seeker  or  guidance  system. 

2.  Missile  homing  to  impact. 

For  an  individual  missile  to  successfully  engage  a  target,  it  must  complete  both 
phases.  Expendable  decoys  can  disrupt  either  of  these  phases.  If  the  missile  cannot  be 
distracted  during  the  target  acquisition  phase  then  it  might  be  "  seduced  "  in  its  homing 
process  (  breaking  the  acquisition  of  the  missile  seeker  ).  The  missile  might  then  be  di¬ 
verted  onto  another  ship,  or  an  empty  zone  (  occupied  by  no  ships),  or  the  diverted 
missile  might "  splash-down  "  to  the  ocean  {  missile  might  divert  not  only  from  its  course 
but  also  it  might  divert  to  an  altitude  closer  to  the  ocean  level  ).  Hence  the  missile  might 
"  splash  "  into  the  ocean  as  a  result  of  losing  its  altitude  after  seduction.  Also,  the  ship 
air  defense  guns  might  shoot  down  the  missile.  If  the  missile  is  not  seduced  then  it  might 
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be  shot  down  by  ship  air  defense  guns,  or  if  it  is  not  shot  down  by  ship  guns  it  may  hit 
a  ship  in  the  BG. 

For  many  missiles  vs  many  targets  (  ships  ),  it  is  assumed  that  the  missiles  act  in¬ 
dependently  and  their  errors  are  not  correlated.  It  is  also  assumed  that  the  effects  of  the 
decoys  are  independent.  These  assumptions  of  independence  may  hold  in  some  opera¬ 
tional  circumstances  but  in  others  they  may  not  hold. 

The  methods  of  defense  against  ASCM  attacks  for  a  BG  that  has  at  least  one  or 
two  high-value  ships  and  various  combat  and  escort  ships  taken  into  account  are: 

1 .  Hard  kill  by  ship  guns  or  antimissile  missiles. 

2.  Soft  kill  by  using  EW  assets  (  splash  down  or  diversion  by  using  decoys ). 

In  this  study,  we  developed  a  simple  model  that  enables  us  to  compare  different 
factors,  such  as  effectiveness  of  different  types  of  decoys  (  seduction  and  distraction  ), 
the  effect  of  changing  the  numbers  and  types  of  decoys,  decoys  vs  hard  kill,  and  non¬ 
linear  interactions  caused  by  diversion  onto  friendly  ships  in  a  multiple-ship  BG. 

In  this  study,  the  following  conditions  hold  for  ASCM  attacks  on  a  multiple-ship 
BG  environment  situation: 

1.  Effectiveness  is  measured  as  survivability  of  the  entire  BG;  not  only  for  a  single 
ship. 

2.  Multiple-ship  environment;  one  or  two-high  value  ships  such  as  aircraft  carriers, 
battle  ships,  or  supply  ships  in  or  near  the  center  of  the  BG.  Low-value  ships,  such 
as  destroyers,  are  far  from  the  center  of  the  BG. 

3.  Xonlinearities;  ship- missile- ship  interactions  are  important,  not  just  one  ship- 
missile/missiles  interaction. 

4.  Decoys  can  be  expendable  decoys  deployed  from  the  ships  for  defending  the  ship 
itself,  or  the  BG,  against  incoming  ASCM  attacks.  Two  types  of  expendable  decoy  s 
are  considered: 

a.  Distraction  decoys  provide  alternate  targets  for  the  missile  during  the  target 
acquisition  process. 

b.  Seduction  decoys  disrupt  the  homing  process. 

Since  this  is  an  unclassified  study,  some  additional  assumptions  regarding  perform¬ 
ance  parameters  are  made  to  exercise  the  model.  The  model,  computer  code,  and  sensi¬ 
tivity  studies  can  be  readjusted  for  special  conditions.  Hence,  the  model  and  sensitivity 
studies  give  results  that  depend  on  unclassified  information. 


2 


II.  BACKGROUND 


A.  BATTLE  GROUP  (  BG  ) 


As  depicted  in  Figure  1,  an  idealized  BG  is  a  group  of  ships  that  has  at  least  one 
or  two  high-value  ships  and  other  low-value  ships. 


. 
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Figure  1.  An  idealized  battle  group 


The  approximate  number  of  total  ships  in  the  group  is  8-15,  which  can  be  divided  into 

two  categories  of  ships  based  on  their  importance.  These  categories  are: 

1.  High-value  ships 

2.  Low-value  ships 
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High  value  ships  may  be  aircraft  carriers  (CV's),  Battle  ships  (BB's),  and  supply 
ships  (AOE's).  Low  value  ships  may  be  guided  missile  cruisers  (CG's),  destroyers  (DD's), 
and  frigates  (FF's).  The  values  of  the  ships  in  a  BG  may  be  related  to  the  mission  of  the 
BG  ,  location  of  the  operation  scene,  and  other  factors  that  affect  the  BG  scenario. 

The  basic  weapons  of  a  BG  for  defending  itself  against  ASCM  attacks  are: 

1.  Guns. 

2.  Antimissile  missiles. 

3.  Electronic  Warfare  assets,  such  as  electronic  counter  measures 

(  ECM  ),  electronic  counter  counter  measures  (ECCM),  and  decoys  for  seduction  or 
distraction. 

In  this  study  we  are  primarily  interested  in  the  weapons  for  hard  kills  (  guns,  anti¬ 
missile  missiles  )  and  decoys  for  soft  kills. 

Mission  types  for  a  BG  usually  require  it  to  remain  on  the  ocean  for  a  long  period 
of  time  and  to  conduct  some  independent  operations.  Generally,  the  main  threats  to  a 
BG  come  from  the  air,  under  the  water,  and  the  surface.  But,  in  modern  warfare,  the 
most  likely  unexpected  method  of  attacking  a  BG  can  be  from  under  water  or  air.  Today- 
most  conventional  submarines  may  use  ASCM's  to  attack  surface  targets.  Hence  we 
discuss  ASCM  attacks  from  different  platforms. 

B.  ANTISHIP  CRUISE  MISSILES  (  ASCM'S  ) 

ASCM  attacks  on  a  BG  may  come  from  any  direction  relative  to  the  course  of  the 
BG.  ASCM's  can  follow  different  flight  patterns.  Two  possible  patterns  can  be  seen  in 
Figure  2.  The  processes  that  an  ASCM  goes  through  in  its  flight  pattern  are: 

1.  Target  acquisition  process. 

2.  Homing  process  to  impact. 

The  target  acquisition  process  for  an  ASCM  includes  tracking,  locating,  identifica¬ 
tion  and  classification  of  the  target;  calculations  for  fire  control;  launching  the  missile; 
flight  until  the  final  phase,  which  includes  precise  tracking  to  the  prelocated  position  of 
the  target  (  homing  to  impact  ).  The  target  acquisition  process  can  vary  for  different 
types  of  missiles.  We  assume  that  this  process  occurs  approximately  from  launching  to 
6000  -  20000  yards  from  the  target  that  the  anti-ship  cruise  missile  (  ASCM  )  is  aimed 
at.  The  range  of  this  process  for  a  missile  depends  upon  the  characteristics  of  the  guid¬ 
ance  system  of  the  missile. 
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Figure  2.  Target  acquisition  and  homing  processes  of  different  ASCM's 


The  homing  process  to  impact  for  an  ASCM  starts  at  the  end  of  the  target  acquisi¬ 
tion  process  and  continues  to  impact  or  miss  of  the  target.  The  homing  process  occurs 
in  relatively  short  distances  compared  to  the  target  acquisition  process.  Although  it 
varies  from  missile  to  missile,  we  assume  that  homing  process  can  occur  around  6000  - 
20000  yards. 

The  speed  of  the  incoming  missile  can  vary  from  mach  1.5  to  mach  2.5  for  different 
types  of  missiles. 

C.  INTERACTIONS,  ASSUMPTIONS,  AND  THE  NATURE  OF  THE  PROBLEM 

The  speed  of  the  BG  is  usually  from  10  to  30  knots.  Considering  the  area  that  is 
covered  by  a  BG,  and  comparing  the  speeds  of  the  BG  and  the  ASCM,  we  can  assume 
that  because  the  speed  factor  of  the  BG  is  small  compared  to  the  speed  factor  of  the 
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incoming  missiles,  speed  does  not  have  a  major  effect  on  avoidance  maneuvers  of  the 
BG  against  ASCM's.  For  this  reason  we  can  assume  that  the  BG  is  almost  fixed  during 
the  attack.  Hence,  we  consider  that  the  distance  traveled  by  the  ASCM  in  the  target 
acquisition  process  is  longer  than  the  distance  traveled  by  the  ASCM  in  the  homing 
process.  Thus,  this  assumption  about  the  fixed  location  of  the  BG  is  even  more  justified 
in  the  homing  process  than  the  target  acquisition  process. 

As  previously  mentioned,  a  BG  that  has  different  kinds  and  values  of  ships  can  be 
attacked  by  ASCM's  launched  by  surface  ,  air  or  underwater  platforms.  In  each  case 
ASCM  's  are  the  faster,  most  unexpected,  precise,  and  preferred  manner  of  attack  against 
a  BG.  Each  different  type  of  missile  usually  shows  almost  the  same  flight  characteristics. 

Hard  kill  devices  of  the  BG  for  defending  itself  against  ASCM's  are  high  speed  small 
caliber  air  defense  guns,  such  as  20  millimeters,  35  millimeters  and  40  millimeters  —  for 
example,  the  PHALANX  system.  Short  range  surface  to  air  missiles  (  SAM's )  are  also 
used  —  for  example  the  SEA  SPARROW  air  defense  missile  system.  Since  this  study 
contains  unclassified  information,  a  hypothetical  hard  kill  probability  (  PHK  )  of  the  BG 
for  an  incoming  missile  is  chosen. 

Soft  kill  devices  (  decoys  )  are  of  two  types:  distraction  decoys  and  seduction  decoys. 
The  quality  factor  (  Q  )  for  distraction  decoys  can  be  defined  as  a  measurement  of  the 
distraction  ability  of  a  distraction  decoy.  If  a  distraction  decoy  has  Q  =  1.0,  it  has  the 
same  effect  as  a  real  ship  on  a  missile  e.g.  essentially  the  same  Radar  Cross  Section  ( 
RCS  ).  Quality  factors  take  on  values  in  the  interval  [0,1] 

The  probability  of  breaking  the  lock  of  a  missile  by  a  seduction  decoy,  (  B  ),  can 
be  defined  as  a  measurement  of  the  seduction  ability  of  a  seduction  decoy.  After  se¬ 
ducing  an  incoming  missile,  if  the  seduction  decoy  successfully  breaks  the  lock  of  the 
incoming  missile,  the  missile  may  divert  or  splash  into  the  ocean.  For  different  flight 
patterns  there  are  different  splash  probabilities.  For  example,  the  flight  pattern  in  Fig. 
2. a  has  a  higher  splash  probability  than  the  one  in  Fig.  2.b 
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III.  SURVIVABILITY  OF  BATTLE  GROUP 


A.  THE  MODEL  TO  CALCULATE  THE  PROBABILITY  MASS  FUNCTION 
VALUES  (  PMF  )  FOR  NUMBER  OF  ENCOUNTERS. 

An  idealized  BG,  has  approximately  8-15  ships.  If  we  define  N0,  A,  and,  R  by 

1.  Ar0 :  number  of  ships  in  the  BG, 

2.  A  :  area  covered  by  the  BG, 

3.  R  :  radius  of  the  BG, 

then  the  ship  density  (p)  for  the  BG  which  covers  area  A  is 


It  seems  reasonable  to  suppose  the  ships  will  maintain  their  distance  from  each  other 
when  an  attack  is  imminent.  Hence  for  simplicity  of  analysis  we  assume  the  ships  are 
located  on  the  intersections  of  the  grid  lines  within  the  area  A.  Later  we  will  allow  the 
ships  to  be  uniformly  located  on  approximately  parallel  lines  normal  to  the  missile  tra- 

A 

jectory.  The  area  per  ship  in  the  BG  is  clearly  ( ~rr- )  which  can  be  visualized  as  a  square 
j  a  A0 

of  side,  (  j-~  )  ,  see  Figure  3. a.  Hence  the  approximate  distance  (  D  )  between  ships 
is 

■ 

As  seen  in  Figure  3.b,  the  random  variable  L  is  the  distance  which  the  missile  travels  as 
it  traverses  the  BG.  The  number  (  N  )  of  possible  ship  encounters  by  the  missile  can  be 
approximated  as  the  integer  value  closest  to  N  where, 
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Lower  case  n,  1,  and,  n'  will  denote  particular  realizations  of  the  r.v.'s  N,  L,  and,  N'  Since 
2R  is  the  maximum  distance  which  a  missile  can  travel  within  the  BG,  the  maximum 
number  of  encounters  will  occur  for  L  =  2R;  this  is  given  by 


We  assume  the  missiles  meet  at  least  one  line  of  ships  in  the  BG.  In  other  words  we  are 
not  concerned  about  missiles  which  miss  the  BG.  Hence  the  range  of  N  is, 

IZNUlRyPjf . 

In  the  absence  of  intelligence  reports  and  any  detection  of  enemy  missiles,  a  missile  at¬ 
tack  may  come  from  any  direction.  Since  there  is  a  circular  symmetry  because  of  random 
direction  of  attacks,  we  assume  the  BG  has  a  circular  shaped  formation.  A  circular  for¬ 
mation  can  give  more  layer  defense  to  the  high-value  ships  in  the  BG  by  locating  them 
in  or  near  the  center  in  case  of  missile  attacks  from  any  direction,  and  it  can  give  efficient 
ultra  high  frequency  (  LHF  )  '  visible  communication  conditions  among  all  ships  in  the 
BG  and  the  BG  commander. 

In  a  circular  formation,  the  area  (  A  )  covered  by  the  BG  is 

A  =  nR2  . 


So  the  maximum  numbor  of  possible  encounters  is 


A' 


(3.1) 


and  for  any  distance  L 


(3.2) 


Let  the  closest  point  of  approach  (  CPA  )  to  the  center  of  the  BG  for  an  ASCM  be 
the  random  variable,  Y,  which  is  uniformly  distributed  in  |  0 ,R  |  From  Figure  4  we  see 
that  Y,  determines  L,  and  therefore  also  N. 

Let  N"  be  a  random  variable  which  is  X  rounded  to  the  nearest  integer.  Then  we 
can  find  its  probability  mass  function  values  (  PMF  )  based  on  the  underlying  uniform 
distribution  of  Y,  for  any  given  A'0  ( the  number  of  ships  in  the  BG  );  see  Figure  4. 
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Figure  4.  Drawing  for  calculation  of  the  eqn.  3.2  and  3.3 


From  Figure  4  we  sec  that 


hence,  using  eqns.  (3.1),  (3.2),  and,  (3.3) 

If  we  take  the  square  of  the  each  side, 

4  OV0K2-Ao>'2) 

A  =  — - - - - 

R 2  71 


)-( 


4A~0K2 
nR 2 


(3-3) 


Solving  for  Y  \vc  find, 
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_  J 


■2 


!*-(■ 


4iVn 


Y=Rj(\-{ 


N2n 

4iVn 


)) 


Replacing  the  r.v.'s  Y  and  N  with  y  and  n  we  obtain  the  correspanding  relations  between 
the  realizations  of  the  r.v.'s.  The  latter  equations  are  used  to  obtain  the  derivative  in  the 
transformation  formula 


hence,  using fY(y)  =—  , 


—Rnn 


4*0./  1 


n2  7i 

lAo- 


and  finally 


A<»)- 


nn 


(3.4) 


Part  of  the  computer  program  winch  we  discuss  in  Chapter  4  computes  the  proba¬ 
bility  mass  values  (  PMF  )  for  a  discrete  number  of  encounters  using  the  PDF  of  n  wrhich 
was  calculated  in  eqn.  (3.4).  As  seen  in  Figure  5,  the  dotted  line  shows  the  points  which 
are  calculated  by  using  eqn.  (  3.4  )  for  A'0  =  14.  Since  there  are  a  discrete  number  of 
layers  in  reality  we  find  the  probability  mass  values  (  PMF  )  for  each  discrete  layer  by 
calculating  the  sub  portions  of  area  under  the  dotted  line  which  are;  0.5  to  1.5  for  1st. 
layer,  1.5  to  2.5  for  2nd.  layer,  2.5  to  3.5  for  3rd.  layer.  For  Figure  5  the  maximum 
number  of  layers  =  4.2  w'hich  is  calculated  using  the  eqn.  (  3.1  )  for  N0  =  14  in  the  BG. 
Since  the  maximum  X'  is  the  integer  number  closest  to  n,  which  is  4,  we  then  use  3.5  to 
4.2  for  the  probability  of  the  4th.  layer. 

Now  let  us  name  the  small  portions  of  the  area  under  ff{n)  as,  AUA2,  Aif...,  At,  wrhere 
i  is  the  closest  integer  value  to  nm„.  If  we  want  to  calculate  the  area  under  the  dotted  line 
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(  PDF  of  n  ),  which  will  be  called  P  in  this  calculation,  we  can  first  calculate  the  A0  which 
is  the  area  under  the  /A(n)  from  0  to  0.5, 


f0.S  f0.5 

4)  =  J0  fd")dn  =  j0 


Then  we  can  find  the  area  of  probability  for  PDF  of  n  between  0.5  and  n, 
ing,  ( since  N  ^  1,  we  eliminate  probability  associated  with  N'  =  0.) 

P=\-A0 


or  we  can  find  total  area  by  calculating 

P  =  Ax  +  A2  +  A3  +  ...  +  Al  , 


where 


At 


i 

by  calculat- 
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NUMBER  OF  SHIPS  IN  THE  BG  (N)  -  14 


NUMBER  OF  LAYERS  OF  SHIPS  IN  THE  BG 


Figure  5.  Calculation  of  the  PMF  values  for  #  encounters  per  missile 


Either  of  these  two  approaches  should  give  the  correct  result  for  total  probability  under 
the  PDF  of  n>  0.5.  We  shall  use  these  two  approaches  to  check  the  results  of  the  com¬ 
puter  program  presented  in  Chapter  4. 

After  finding  the  total  area  under  the  PDF  of  n^  0.5,  we  normalize  probability  mass 
values  of  each  layer  with; 


PMF,=^- 


PMF^  — 
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PMF, 


and  for  ith.  interval, 

PMF,  =  . 

Now  the  sum  of  all  the  probabilities  should  be 

PMF{  +  PXIF2  +  PMF2  +  ...  +  PMFt  -  1  , 

so  that  the  number  of  possible  encounters  is  an  integer  ^  1  rather  than  >  0. 

Figure  5  gives  an  example  calculation  of  the  normalized  PMF  values  for  the  number 
of  possible  encounters  when  N0  =  14.  The  dotted  line  shows  the  values  of f^n)  and  the 
solid  line  connects  the  normalized  PMF  value  points  for  each  discrete  layer. 


* 
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B.  THE  MODEL  TO  CALCULATE  THE  PROBABILITIES  OF  SHIP  HIT, 
MISSILE  KILL,  MISSILE  DIVERSION,  AND  SHIP  SURVIVABILITY 
1.  Probability  calculation  for  one  layer 

The  probabilities  of  ship  hit ,  missile  kill,  missile  diversion,  and  ship  survivability 
can  be  calculated  by  constructing  a  sample  probability  tree  model  for  one  layer.  As  seen 
in  the  tree  diagram  in  Figure  6,  when  a  missile  arrives  to  a  line  of  possible  encounters  it 
may  meet  one  or  more  ships  or  no  ships.  We  now  assume  that  the  ships  are  located 
randomly  along  the  line  perpendicular  to  the  missile  trajectory. 


dMn/  *««*  NesiHL,  D!v,rtloa 

X  . - .  »  UK 

!  \  TJP  / 

/  'sdMradicL  /  afff  sriaj.fr  - 

at  least  \  /hard  kills 

che  ship  \ Ruction  ^ot£Fiosh 

/  }  \  • — - •  Diversion 

— 1 hal’d  Jdll^  «»»• 

V  seflietlcB  ^ 

X  net  ^  Ship 

■hflrCl  destroyed 

—hal'd  kill ^ 


“L _ < 

arrives  X 
&t  &  line 


nor  hara^rH*1 


Diversion 


Figure  6.  Probability  tree  diagram  for  one  layer  of  encounters 


As  can  be  seen  in  Figure  7,  if  a  is  the  angular  search  beam  of  the  missile,  then  the 
probability  a  missile  does  not  detect  a  ship  on  the  layer  (  PONS  )  is 

-23 


PONS = e 


(3.5) 


Figure  7.  Calculation  for  probability  if  missile  meets  at  least  a  ship 

which  is  based  on  the  Poisson  Approximation  to  the  Binomial  Distribution  of  the  num¬ 
ber  of  ships  falling  within  the  angular  search  beam  width  of  the  missile.  On  the  other 
hand  the  probability  a  missile  detects  at  least  one  ship  per  layer  (  POS  )  is  obviously 

POS  =  1  -  e~u  .  (3.6) 

a.  Missile  meets  at  least  one  ship 

If  a  missile  meets  at  least  one  ship,  then  it  can  be  successfully  distracted  or 
not  distracted  by  the  ship  using  a  distraction  decoy  with  quality  factor  Q.  The  dis- 


traction  probability  (  POD  ),  which  is  the  ratio  of  the  equivalent  number  of  false  targets 
(  QD  )  to  the  total  number  of  targets  (  1  +  QD  )i  ,  and  is  assumed  to  be  given  by: 


POD  =  ( 


QD 


(1  +  QD) 


) 


(3.7) 


where  D  is  the  number  of  distraction  decoys  deployed  by  a  given  ship  during  the  missile 
attack.  On  the  other  hand,  the  probability  of  no  distraction  of  a  missile  given  it  meets  a 
ship  (  POND  )  is 


POND  =  ( 


1 


(1  +QD) 


)  • 


(3.8) 


Then: 

If  distraction  occurs,  then  the  missile  can  be  destroyed  by  a  hard  kill;  The 
probability  of  hard  kill  given  the  missile  is  distracted  is 


PMK,  =  (1  -  e-2"X  )PHK  . 


(3.9) 


If  the  missile  is  not  hard  killed  it  can  splash  down  and  be  destroyed  (  soft  kill )  or  not 
splash  down  and  instead  be  diverted  after  the  distraction  of  the  missile  occurs.  The 
probability  of  splash  down  of  a  missile  =  PSPL.  If  splash  occurs  obviously  we  can  ac¬ 
cept  this  as  a  missile  soft  kill.  On  the  other  hand,  even  if  it  is  not  splashed  it  is  still  di¬ 
verted.  In  this  case  the  probability  of  missile  kill  given  the  missile  is  distracted  is 

PMK2  =  (1  -  e~2g)(  ■  )(1  -  PHK)PSPL  ,  (3.10) 

and  the  probability  of  diversion  (  PDIV  )  given  the  missile  is  distracted  is 

PD1VX  =  (1  -  e~2a)(  )0  -  PHK){\  -  PSPL)  .  (3.11) 

If  a  missile  is  not  distracted  with  the  probability  calculated  in  eqn.  (3.7)  it 
can  be  seduced  or  not  seduced  by  the  BG  using  seduction  decoys.  In  this  case  the 
probability  of  seduction  (  PSED  )  is 
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(3-12) 


PSED  =  (1  —  (1  -  B)M5)  , 
and,  the  probability  of  no  seduction  is 

PONS  =  (1  -  B)ms  (3.13) 

respectively,  where 

1.  B  *  Probability  of  breaking  the  lock  of  a  missile  for  a  seduction  decoy  in  the  BG. 

2.  MS  =  Number  of  seduction  decoys  in  per  ship. 

If  seduction  occurs  the  missile  can  be  destroyed  by  hard  kill  or  not  destroyed  by  hard 
kill.  The  probability  of  missile  is  hard  killed  given  the  seduction  occurs  is 

PMK3  =  (1  -  e~23)(  J^QD)  )(1  "  (1  “  B)MS)PHK  .  (3.14) 

The  probability  of  no  hard  kill  given  that  seduction  occurs  is 

PMK,  =  (1  -  e~u)(  )(1  -  (1  -  Bfs)(  1  -  PHKjPSPL  ,  (3.15) 

and  the  probabability  of  diversion  given  the  seduction  occurs  is 

PDIV2  =  (1  -  e~2*)(  -  +lQD)  )(1  -  (1  -  g)"5)(l  ~  PSPL)  .  (3.16) 

If  seduction  does  not  occur,  then  the  probability  of  hard  kill  given  the  seduction  does 
not  occur  is 


PMK5  =  (1  -  e~2a)(  {l'+QD)  X0  -  B)ms)PHK  ,  (3.17) 

or  if  it  cannot  be  hard  killed  it  may  hit  a  ship  so  the  probability  of  a  missile  hitting  a  ship 
given  the  seduction  does  not  occur  is 

PSHI  =  (1  -  e~2a)(  )((1  -  5)M5)(1  -  PHK)  .  (3.18) 


b.  Missile  meets  no  ship 


If  a  missile  meets  no  ships,  then  it  can  still  be  destroyed  by  hard  kill,  then 
the  probability  of  hard  kill  given  the  missile  meets  no  ships  is 

PMK6  =  e~2aPHK  ,  (3.19) 

and  the  probability  of  diversion  given  the  missile  meets  no  ships  is 

PDIV3  =  e~2a(l  -  PHK)  .  (3.20) 

From  all  of  these  equations  we  calculated  above,  we  can  calculate  the  total 
missile  kill  probability  (  TPMK  )  per  one  line  of  encounters  by  using  the  equations. 

TPMK  =  PMKy  +  PMK2  +  PMK3  +  PMKA  +  PMKS  4-  PMKt  (3.21) 

If  we  write  the  whole  equation  for  TPMK; 

TPMK  =  ( 1  -  e~2x)(  — )PHK  +  ( 1  -  e~2*)(  )( 1  -  PHK)PSPL 

+  (1  -  ")( TiTm  )(l " (1  -  B)MS)PHK + (1 '  ^ UTm } 

(1  -  0  -  B)MS){  1  -  PHK)PSPL  +  (1  -  e~2a){  - }  )((1  -  B)"S)PHK 

+  {e~2x)PHK  .  (3.22) 

The  probability  of  ship  hit  in  the  BG  (  PSHI  )  is  again, 

PSHI  —  (1  -  g~2a)(~(i ~+QD)m  "  B)MS){{  ~  PHK)  ’  (3  23) 

so  survivability  of  BG  (  the  probability  that  all  ships  survive  the  attack  )  as  a  measure 
of  effectiveness  for  one  line  (  PSSUR  )  is 

PSSUR  =  (I  —  (1  —  e_2a)(  (1  +1qD)  )(1  -  5)WS(I  -  PHK))  .  (3.24) 

The  total  diversion  probability  (  TPDIV  )  is  then 
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TPDIV=  PDIV]  +  PDIV2  +  PDIV2 


or 


TPDIV  -  (1  -  e~2a)(  -(]  fDQ--  )(1  -  PHK){\  -  PSPL )  +  (1  -  e~2*){  (i  +1-Q/))  ) 

(1  -  (1  -  5)W5)(1  -  PSPL)  +  ( e~2*)PHK  .  (3.25) 

Equations  3.21,  3.22,  3.23,  3.24  are  calculated  on  the  probability  tree  dia¬ 
gram  for  one  line  of  encounters  in  Appendix  B.  The  reader  can  refer  to  Appendix  B  to 
get  an  overall  idea  for  understanding  better. 


2.  Probability  calculation  for  ail  the  layers 

We  have  shown  how  to  calculate  the  maximum  number  of  lines  using  eqn.  (3.2), 
the  probabilities  of  any  ship  being  hit  in  the  BG,  missile  kill,  missile  divert,  and 
survivability,  for  one  line  of  encounters.  Now  we  can  find  the  probabilities  conditioned 
on  the  number  of  layers  (  or  potential  encounters  with  ships )  which  make  up  the  whole 
BG.  For  these  calculations  we  shall  use  a  probability  tree,  which  can  be  seen  in  Figure 
8,  for  all  the  layers  in  the  BG,  and  we  will  be  using  the  equations  that  we  calculated  for 
one  layer  encounters. 

As  we  can  see  in  Figure  7,  for  each  line  the  missile  may  hit  a  ship  or  get  diverted 
ot  get  "  killed  ".  If  a  missile  hits  a  ship  or  gets  killed,  the  threat  of  the  missile  is  termi¬ 
nated.  But  if  a  diversion  occurs  the  same  conditions  and  probability  calculations  take 
place  again  and  again  until  it  passes  the  last  line  of  encounters  within  the  BG.  For  the 
last  line,  diversion  is  equivalent  to  the  termination  of  the  missile. 

By  using  the  equations  in  Part  1  and  also  using  Figure  7  and  Figure  8  we  can 
find  the  following  probabilities  for  a  BG: 

Probability  of  missile  destruction  for  a  BG  (  PMKBG  )  averaged  over  the 
number  of  lines  of  encounters  for  a  given  number  (  N  )  of  ships  is 


PMKBG  =  YjPr  iMK  |  i  layers]  PMF,  .  (3.26) 

/=  i 

For  instance  by  adding  the  probabilities  for  mutually  exclusive  (  M.E.)  events  we  have: 
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PrlMK  |  3  layers]  =  PMK  +  (PDIl)  PMK  4-  {PDIV7)  PMK  . 

Probability  of  a  ship  getting  hit  for  a  BG  (  PSUBG  )  which  has  n  lines  of  en¬ 
counters  for  a  given  number  (  N  )  of  ships  is 


PSHBG  =  YjP>  iSI{I  I '  /«)*«]  pMpt  •  (3.27) 

/=i 

For  instance  by  adding  probabilities  for  M.E.  events  we  have, 

Prism  |  3  layers']  =  PS  11 1  +  (PDHr)  PS  111  4-  {PDIV7)  PS  II I  . 
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Probability  of  diversion  of  a  missile  for  a  BG  (  PD1VBG  )  which  has  n  lines  of 
encounters  for  given  number  (  N  )  of  ships  is 

n 

PDIVBG  =  Y^PDIV)1  PMF ,  .  (3.28) 

i=  1 

Also  we  can  calculate  the  probability  of  ship  survivability  (  PSSBG  )  as  measure 
of  effectiveness  (  MOE  )  for  the  BG. 

PSSBG  =  (1  -  PSHBG)  .  (3.29) 

We  will  be  using  the  probabilities  we  calculate  in  eqns.  3.26,  3.27,  3.28  and  3.29 
for  drawing  the  graphs  to  understand  the  different  relationships  between  the  survivability 
of  the  BG  and  factors  such  as  the  number  of  each  different  type  of  decoy,  the  number 
of  ships  in  the  BG,  the  number  of  missiles  attacking  the  BG,  the  distraction  decoy  Q 
factor,  the  seduction  decoy  B,  factor,  etc. 

To  illustrate  these  relationships  we  can  calculate  probabilities  in  eqns  .3.26,  3.27, 
3.28  and  3.29  by  using  the  FORTRAN  77  program  in  the  Appendix.  After  we  run  the 
program  by  changing  the  different  parameters  we  get  the  example  results  as  Figures  9  - 
17. 

As  seen  in  Figure  9  we  run  the  program  with  the  parameters: 

1.  Number  of  ships  in  the  BG  (A'0)  =  5. ..50. 

2.  Number  of  distraction  decoys  per  ship  (  DI  )  =  1. 

3.  Number  of  seduction  decoys  per  ship  (  MS  )  =  1. 

4.  Hard  kill  probability  of  BG  (  PHK  )  =  0.25. 

5.  B  factor  of  seduction  decoys  =  0.2. 

6.  Q  factor  of  distraction  decoys  =  0.2. 

7.  Splash  probability  of  missile  (  PSPL  )  =  0.05. 

8.  Alpha  angle  of  missile  (a)  =  15’. 

and  later  we  get  the  ship  survival  probabilities  as  the  number  of  ships  in  the  BG  are  in¬ 
creasing.  We  can  easily  see  from  the  graph  that  when  the  number  of  ships  increases  in 
the  BG,  survivability  of  ships  decreases.  Also  ship  survivability  of  BG  has  a  more  nega¬ 
tive  slope  {  decreasing  survivability  )  for  N0  -  5  to  15  than  N0  =  15  to  50  (  50  can  be 
thought  of  as  very  large  =  infinite  number  of  ships  for  BG  ).  When  the  number  of  ships 
goes  to  infinity,  then  the  probability  of  survivability  of  BG  is  almost  constant. 
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SHIP  SURVIVABILITY  IN  THE  BG 
N~5~50,  DI-MS- 1,  Q-B-.2,  PHK-.25 
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Figure  9.  Probability  of  no  hits,  while  number  of  ships  varies 


If  we  look  at  Figure  10  and  the  Figure  11  we  can  easily  see  that  the  different 
types  of  decoys  affect  the  survival  probability  differently.  In  Figure  10: 

1.  The  number  of  ships  (AT0)  varies  from  5-30. 

2.  Number  of  distraction  decoys  (  DI )  varies  from  1-5. 

We  get  five  different  survival  curves  for  BG.  while  the  other  parameters  remained  un¬ 
changed  and  are, 

3.  Number  of  seduction  decoys  per  ship  (  MS  )  =  1. 

4.  Probability  of  hard  kill  level  of  BG  (  PHK  )  =  0.25 

5.  Splash  probability  of  the  missile  (  PSPL  )  =  0.05 

6.  B  factor  of  seduction  decoys  =  0.2 
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SURVIVABILITY  OF  BG 

N^5-30,  DI-1-5,  MS-1,  PHK-.25,  PSPL-.05,  B-Q-.2 


Figure  10.  Probability  of  no  hits,  distraction  decoys  vary  from  1  to  5 

7.  Q  factor  of  distraction  decoys  =  0.2 

8.  Alpha  angle  of  the  missile  (a)  =  15' 

In  Figure  11: 

1.  Number  of  distraction  decoys  (  DI )  =  1. 

2.  The  number  of  seduction  decoys  (  MS  )  varies  from  1-5. 

We  get  five  different  survival  curves  for  the  BG.  while  the  other  parameters  remained 
unchanged  and  are  the  same  as  above.  As  we  see  in  Figures  10  and  11,  different  types 
of  decoys  affect  the  survivability.  In  this  case  the  number  of  seduction  decoys  has 
slightly  more  effect  than  the  distraction  decoys  on  survivability  of  a  BG. 
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SURVIVABILITY  OF  BG 


N-5-30,  MS-1 -5,  01-1,  PHK-.25,  PSPL-.05,  B-Q-.2 


Figure  1 1.  Probability  of  no  hits,  seduction  decoys  vary  from  1  to  5 

If  we  change  the  number  of  decoys,  but  we  do  not  change  the  number  of  ships 
in  the  BG,  we  can  get  Figure  12  and  Figure  13.  In  Figure  12: 

1.  Number  of  distraction  decoys  (  DI  )  varies  from  1-5. 

2.  The  number  of  ships  (Nt)  =  5. 

3.  Number  of  seduction  decoys  per  ship  (  MS  )  =  1. 

4.  Probability  of  hard  kill  level  of  BG  (  PHK  )  =  0.25. 

5.  Splash  probability  of  the  missile  (  PSPL  )  =  0.05. 

6.  B  factor  of  seduction  decoys  =  0.2. 

7.  Q  factor  of  distraction  decoys  =  0.2. 

8.  Alpha  angle  of  the  missile  (a)  =  15'. 
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Figure  12.  Survivability  of  BG  while  num.of  distraction  decoys  varies 


In  Figure  13,  we  modify  the  first  and  third  parameters  so  that: 

1.  Number  of  distraction  decoys  (  DI  )  =  1. 

3.  Number  of  seduction  decoys  (  MS  )  varies  from  1  to  5. 

In  these  two  figures,  when  the  seduction  decoys  were  changed,  we  obtain  a  slightly 
higher  survival  probability  compared  to  when  the  number  of  distraction  decoys  were 
changed.  The  benefit  of  expendables  are  apparent  from  these  figures. 

The  number  of  missiles  attacking  the  BG  is  also  an  important  factor  for 
survivability,  which  is: 

Pr  (  Survival  given  M  missile  attacks  )  =  Pr{  Survivability  of  BG  )M  ,  M  =  1,2,... 
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SURVIVABILITY  OF  BG 

N-15,  MS- 1-5,  01-1,  PHK-.25,PSPL-.05,B-Q«.2 


Figure  13,  Survivability  of  BG  while  num.of  seduction  decoys  varies 

This  assumes  independence  for  each  missile  and  that  the  decoys  are  not  depleted  during 
successive  attacks.  After  running  the  program  and  calculating  the  survival  probabilities 
for  M  missiles  we  got  the  graphs  in  Figure  14  and  Figure  15.  In  Figure  14  we  run  the 
program  for  M  missile  attacks  while, 

1.  Number  of  missile  attacking  the  BG  (  M  )  varies  from  1  to  10. 

2.  Number  of  distraction  decoys  (  DI  )  varies  1  to  5. 

3.  Number  of  seduction  decoys  (  MS  )  =  I. 

4.  The  number  of  ships  (AT0)  =  10. 

5.  Probability  of  hard  kill  level  of  BG  (  PHK  )  =  0.70. 

6.  Splash  probability  of  the  missile  (  PSPL  )  =  0.05. 
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SURVIVABILITY  OF  BG  IN  M-1-N  MISSILE  ATTACKS 
Nf  10,  01-1-5,  MS-1,  B-.5,  Q— .7,  PHK-.7,  PSPL-.05 


Figure  14.  Survivability  of  BG  while  num.of  distraction  decoys  varies 

7.  B  factor  of  seduction  decoys  =  0.5. 

8.  Q  factor  of  distraction  decoys  =  0.7. 

9.  Alpha  angle  of  the  missile  (a)  =  15”. 

In  Figure  15  we  run  the  program  for  M  missile  attacks  making  the  following  changes: 

1.  Number  of  seduction  decoys  (  DI  )  varies  I  to  5. 

2.  Number  of  distraction  decoys  (  MS  )  =  1. 

From  Figure  15  and  16,  we  can  see  that  seduction  decoys  have  more  effect  on 
survivability  than  distraction  decoys  have.  Especially  after  the  third  decoy  has  been  de¬ 
ployed.  For  each  type  of  decoy  there  is  not  a  very  big  difference  on  survivability  for 
number  of  attacking  missilies  (  M  )  =  1,  but  when  M  increases  the  seduction  decoys 
have  more  effect  yielding  significantly  higher  survivability. 

The  specific  values  of  Q  and  B  are  important  factors  for  determining  the  effect 
of  expendibles.  Figures  16  and  17  show  the  effects  of  varying  the  Q  factor  of  a  dis- 
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Figure  15.  Survivability  of  BG  while  num.of  seduction  decoys  varies 

traction  decoy  and  the  B  factor  of  a  seduction  decoy,  for  the  case:  1  seduction  and  1 
distraction  expended  for  engagement  by  1  missile.  In  each  case  B  and  Q  varied  between 
0.0  to  1.0  while  the  other  parameters  remained  unchanged.  The  effect  of  seduction  decoy 
B  factor  to  survivability  is  higher  than  the  effect  of  distraction  decoy  Q  factor.  In  surface 
ships  Q  and  B  may  assume  any  value  between  0  to  1.  The  specific  values  of  B  and  Q 
depend  on  systems  employed  and  the  dynamics  of  the  engagement.  In  Figure  16,  we  run 
the  program  for  different  values  of  distraction  decoy  Q  factor  while  other  parameters 
remain  unchanged. 

1.  Q  factor  of  distraction  decoys  varies  from  0  to  1. 

2.  Number  of  seduction  decoys  (  MS  )  =  1. 

3.  The  number  of  ships  (Ar0)  =  15. 

4.  Number  of  distraction  decoys  (  DI  )  =  1. 

5.  Probability  of  hard  kill  level  of  BG  {  PHK  )  =  0.25. 

6.  Splash  probability  of  the  missile  (  PSPL  )  =  0.05. 
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SURVIVABILITY  OF  BG 

N— 10,  Q-0--1.  B-.2,  DI-MS-1,  PHK-.25,  PSPL-.05 


DISTRACTION  DECOY  Q  FACTOR 


Figure  16.  Survivability  of  BG  while  Q  factor  of  distraction  decoys  varies 

7.  B  factor  of  seduction  decoys  =  0.2. 

8.  Alpha  angle  of  the  missile  (a)  =  15'. 

In  Figure  17  we  run  the  program  for  different  values  of  the  seduction  decoy  B  factor 
while  other  parameters  stay  unchanged. 

1.  B  factor  of  seduction  decoys  varies  from  0  to  1  (  0.1  increments ). 

2.  Q  factor  of  distraction  decoys  =  0.2. 
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SURVIVAL  PROBABILITIES 


SURVIVABILITY  OF  BG 

R-10,  B-O-vl.O  ,Q-.2,  MS-DI-1,  PHK-.25,  PSPL-.05 


SEDUCTION  OECOV  B  FACTOR 


Figure  17.  Survivability  of  BG  while  B  factor  of  seduction  decoys  varies 


31 


IV.  CONCLUSION 

While  defending  the  BG  against  anti-ship  cruise  missile  attacks  in  the  open  ocean, 
a  major  concern  is  survivability  of  the  BG  and  the  impact  of  the  different  performance 
factors  on  the  survivability  of  the  BG. 

Since  there  is  a  big  difference  between  the  speed  factors  of  a  BG  and  a  missile,  we 
have  assumed  that  the  BG's  location  is  fixed  in  this  particular  situation.  We  also  defined 
the  measure  of  effectiveness  (  MOE  )  as  the  probability  that  all  ships  survive  the  attack 
of  M  missiles. 

Major  factors  which  impact  the  BG's  survivability  are  number  of  ships  in  the  BG, 
the  BG's  hard  kill  level,  and  the  BG's  expendable  decoys.  W'e  have  seen  the  benefit  of 
expendables  in  the  examples  of  the  previous  chapter  which  were  drawn  using  the  data 
calculated  by  the  FORTRAN  77  computer  program  in  the  Appendix.  Specifically,  we 
have  investigated  the  sensitivity  to  the  distraction  decoys  Q  factor  and  the  seduction 
decoys  B  factor  which  are  important  for  determining  the  effects  of  expendables. 

The  analytic  results  of  the  model  which  were  calculated  by  the  computer  program 
give  results  consistent  with  those  expected  in  the  beginning.  If  the  number  of  ships  in¬ 
creases  and  the  other  factors  remain  constant,  survivability  of  BG  /  ship  decreases,  in 
part  because  of  the  existence  of  the  nonlinear  interactions  between  ships  in  the  BG.  The 
tactics  of  single  ship  defense  against  ASCM's,  the  position  of  the  ship  in  the  BG  (  the 
ship's  nearness  to  center  /  high-value  ships  )  may  also  affect  the  BG's  multiple  ship 
overall  survivability,  or  the  survivability  of  high  value  targets  such  as  CV's. 

Different  types  of  decoys  affect  survivability  differently.  This  is  done  by  either  the 
B  or  Q  values  of  the  different  type  of  decoys,  and  /  or  the  number  of  decoys  that  can 
be  deployed.  Specifically  distraction  decoys'  Q  factor  and  seduction  decoys'  B  factor  are 
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important  for  determining  the  effects  of  expendables.  Hence  one  should  analyze 
survivability  of  the  BG  and  include  cost  factors  of  the  expendables  in  the  model  before 
allocating  and  deploying  them  on  each  ship.  Different  trade-offs  also  can  be  made  be¬ 
tween  decoys,  hard  kill  levels,  and  number  /  position  of  ships  in  the  BG,  which  has  lim¬ 
ited  space  and  budget  for  devices,  weapons,  amount  of  decoys,  and  missiles. 

When  we  include  the  cost  factors  of  these  variables  in  our  model  we  can  easily  get 
trade-offs  by  using  a  survivability  cost  comparision  between  different  types  of  assets. 
By  using  the  results  of  this  model  we  can  change  the  different  variables  such  as  the 
number  of  decoys  and  hard  kill  levels,  subject  to  given  cost  constraints,  to  optimize  the 
survivability. 


33 


APPENDIX  A.  FORTRAN  PROGRAM  TO  CALCULATE  THE  SHIP 

SURVIVABILITY  IN  BG 
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PROGRAM  SENGEL 

AUTHOR.  CENGIZ  SENGEL 

NPGS  OPERATIONS  RESEARCH. (  360  ) 


21, AUG, 1988  * 
MOD  7  * 

* 


.»  ■■  ■>  I  t  »  .1.  t  !„,l. 

Tr7CTC7tTCTCTf7CTCyrTC' 


THIS  FORTRAN  77  PROGRAM  CALCULATES  THE  SHIP  SURVIVAL  PROBABILITY 
IN  THE  BATTLE  GROUP  (  BG  )  FOR  N  GIVEN  NUMBER  OF  SHIPS  AND  OTHER 
PARAMETERS  WRITTEN  BELOW.  PROGRAM  ALSO  CALCULATES  THE  SHIP 
HIT  FOR  BG,  SINGLE  SHIP  HIT  /  SURVIVAL,  ATTACKING  MISSILE, 

DIVERSION,  KILL  PROBABILITIES  FOR  M  MISSILE  ATTACKS.  ASSUMING 
INDEPENDENT  MISSILE  ATTACKS  M  CAN  BE  1,2,... ,N  AND  BG  HAS  ALMOST 
CIRCULAR  FORMATION  AGAINST  MISSILE  THREAT  WHICH  MAY  COME  FROM 
ANY  DIRECTION. 

*********** 

PARAMETERS: 

N=  #  OF  SHIPS/INPUT 
L=MAXIMUM  LAYERS  FOR  GIVEN  N  (REAL) 

MXL=  L  VALUE  ROUND  TO  NEAREST  INTEGER 
TSUM=  TOTAL  PROBABILITIES  BETWEEN  . 5  AND  L 
TS  =  TOTAL  PROBABILITIES  WITH  USING  1-PROBABILITY  BETWEEN  0 
AND  .5 

TSUM=  TOTAL  PROBABILITY  AREA  UNDER  PDF  FCN  BETWEEN  0.5  ,MXL 
A=ALPHA  ANGLE  OF  ONCOMING  MISSILE/INPUT 
Q=  QUALITY  FACTOR  OF  DISTRACTION  DECOY  (Q.  LE.  1)/INPUT 
D=  #  OF  DISTRACTION  DECOYS/ INPUT 

B=  PROBABILITY  OF  BREAKING  THE  LOCK  FOR  A  SEDUC.  DECOY/ INPUT 
MS=  #  OF  SEDUCTION  DECOYS/ INPUT 
PHK=  HARD  KILL  PROBABILITY  OF  SHIP  FOR  KILLING  MISSILE/ INPUT* 
PSPL=PROBABILITY  OF  MISSILE  GET  SPLASHED  AFTER  DECOVED/INPUT  * 

PSH=  PROBABILITY  OF  SHIP  GET  HIT  * 

PDV=  PROBABILITY  OF  MISSILE  GET  DIVERTED  * 

PMK=  PROBABILITY  OF  MISSILE  GET  KILLED  (  HARD  KILL  OR  SPLASH)* 

CPSS=  CONDITIONAL  PROBABILITY  OF  SURVIVAL  FOR  BG  (  PROB.  OF  * 

NO  HIT) 

CPDV=  CONDITIONAL  PROBABILITY  OF  MISSILE  DIVERSION 
CPMK=  CONDITIONAL  PROBABILITY  OF  MISSILE  KILL 
CPSS1=  CONDITIONAL  PROBABILITY  OF  SINGLE  SHIP  SURVIVABILITY 
IN  THE  BG 

CPSH=  CONDITIONAL  PROBABILITY  OF  SHIP  HIT  IN  THE  BG 
CPSH1=  CONDITIONAL  PROBABILITY  OF  SINGLE  SHIP  HIT  PROBABILITY* 
IN  THE  BG  * 

CPSSM=  CONDITIONAL  PROBABILITY  OF  SHIP  SURVIVABILITY  IN  CASE  * 

OF  M  MISSILE  ATTACKS  * 

•T7T/1  inTjTTrTrTrir7rjriTinnT#rmT  * 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 
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REAL  N,L,MXL,P ,PR01 ,PL,MXL1 ,P0 ,PI ,TSUM,TS ,A,Q,D,B ,MS ,PSPL,PHK,PSHI 
5 , PSH , PDV , M , PDVI , PMKI , PMK , DI , MXL2 , PL1 , PLO , AL , CPSS , CPSH , CPSS 1 , CPSH1 , 


6CPSS1M,CPSH1M,K,R 

CALL  EXCMS  ( 'FILEDEF  01  DISK  DEMET  OUTPUT  Al') 

DATA  N/10/.A/15.  /.PSPL/.  05/,PHK/. 7/.DI/1. /,Q/. 7/,B /. 5/ 

*  DEFINING  THE  CONSTANT  PI  AND  ALPHA  RADIAN 

PI=3. 141592654 
AL=-2*( A*(PI/180. )) 

DO  14  R=l.  ,5. 

*  DO  14  N=5. ,  50. 

M=N 

MS=R 

*  WRITE(01,24)  N 

24  FORMATC IX, F5. 1) 

♦CALCULATING  THE  MAXIMUM  NUMBER  OF  LAYERS 
Ir=2*SQRT(N/PI) 

♦ROUNDING  THE  MAXIMUM  NUMBER  OF  LAYERS  TO  THE  NEAREST  INTEGER 
MXL=ANINT(L) 

SUM=0 

MXL1=(MXL-1)+.  5 

♦CALCULATION  OF  THE  AREA  UNDER  PDF=TS  (.5,L)  AND  AREA  OF  LAST  PORTION 
P0=ABS( 1-SQRT( l-( .  5**2*PI/(4*N) ) ) ) 

PL=( SQRT( 1.  -(L**2.*PI/(4.  ♦N))))-(SQRT(1. -(MXL1**2. *PI/(4. *N) ) ) ) 

TS=1-P0 

PLP=ABS(PL)/TS 

*  WRITE(01,22)  PLP,L,TS 

22  FORMAT (  IX, F9.  8, 1X.F9.  6,1X,F9. 8) 

T=Q.  0 
CPSH=0. 0 
CPMK=0. 0 
CPDV=0. 0 
PSH=0. 0 
PMKI=0. 0 
PDV=0. 0 

♦CALCULATING  THE  PROBABILITIES  FOR  ONE  LAYER  YOU  CAN  FIND  PMK  BY  TWO  WAY 
PSHI=( l-EXP(AL) )*( ( l/( 1+CQ*DI)))*( ( 1-B)**MS)*( 1-PHK)) 

PDIV=( 1-EXP( AL) )*( ( (Q*DI)/( 1+(Q*DI) ) )*( 1-PHK)*( 1-PSPL) 
1+((1-(1-B)**MS)*(1-PHK)*(1-PSPL)*(1/(1+(Q*DI)))))+ 

2((EXP(AL))*( 1-PHK)) 

♦FIRST  WAY  OF  FINDING  PMK 

*  PMK=l-( PDIV+PSHI) 

♦OR  SECOND  (LONGER  WAY)  OF  FINDING  PMK 

PMK1=((1-EXP(AL))*((Q*DI)/(1+(Q*DI)))*(PHK)) 

PMK2=( ( 1-EXP( AL) )*( (Q*DI ) /( l+( Q*DI ) ) )*( 1-PHK)*PSPL) 

PMK3=( 1-EXP( AL) )♦( l/( 1+(Q*DI ) ) )*( l-( 1-B)**MS)*PHK 

PMK4=( ( 1-EXP( AL) )*( l/( l+( Q*DI ) ) )*( l-( 1-B)**MS)*( 1-PHK)*PSPL) 

PMK5=( ( 1-EXP( AL) )*( l/( 1+(Q*DI) ) )*( ( 1-B)**MS)*PHK)+EXP(AL)*PHK 
PMK6=EXP( AL ) *PHK 
PMK=PMK 1 +PMK2 +PMK3+PMK4+PMK5 
♦CALCULATION  OF  PROBABILITIES  FOR  LAYERS 
DO  15  F=1 ,MXL, 1 
PR=0. 0 

PSH=( PSH+PDIV**(T)*PSHI ) 

T=T+1 

PDV=PDIV**F 
IF(F. EQ. MXL)THEN 
CPSH=CPSH+PSH*PLP 
CPMK=CPMK+PMKI*PLP 
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CPDV=CPDV+PDV*PLP 

ELSE 

PR=SQRT( 1 - ( ( F+.  5 ) **2*P I / ( 4*N ) ) ) - SQRT( 1  ■ -  ( ( F- . 5 ) **2*PI / ( 4*N ) ) ) 

PR01=ABS(PR)/TS 

PMKI=( 1-(PDV+PSH) ) 

CPSH=CPSH+PSH*PR01 
CPMK=CPMK+PMKI*PR01 
C  PDV=C  PDV+PD V*PRO 1 
END  IF 

*  WRITE(01,27)DI,MS 

27  FORMAK  1X,2(F9.  6,  IX)) 

15  CONTINUE 

CPSS=1-CPSH 

♦DISPLAY  CONDITIONAL  PROBABILITY  OF  SHIP  HIT,  MISSILE  KILL,  DIVERSION 

*  WRITE(01,23)CPSH,CPMK,CPDV 
♦DISPLAY  THE  SHIP  SURVIVABILITY  AS  'MOE' 

♦SINGLE  SHIP  PROBABILITY  CALCULATIONS 

*  CPMK1=CPMK/N 

*  CPSH1=CPSH/N 

*  CPSS1=CPSS/N 

♦BG  PROBABILITY  CALCULATIONS  GIVEN  M  MISSILE  ATTACKS  TO  THE  BG 
CPSHM=1.  -(I.  -(CPSH)**M) 

CPSSM=1.  -(1.  -(CPSS)**M) 

♦SINGLE  SHIP  PROBABILITY  CALCULATIONS  GIVEN  M  MISSILE  ATTACKS  TO  THE  BG 
♦WHICH  IS  M  =  1,  2,  3,  ....  N. 

CPSH1M=1. -(1. -(CPSH/M)**M) 

CPSS1M=1. -(1. -(CPSS/M)**M) 

WRITE(01 ,23)N,MS ,DI ,CPSS ,CPSSM 
23  FORMAT (  IX, F4.  1, IX, F4.  1 , 1X,F4.  1 , 1X,2(F9. 7 , IX)) 

14  CONTINUE 

31  STOP 

END 
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APPENDIX  B.  PROBABILITY  CALCULATIONS  FOR  ONE  LINE  OF 

ENCOUNTERS. 

We  can  obtain  all  the  equations  calculated  for  one  line  of  encounters  with  the  aid 
of  the  probabilities  on  the  probability  tree  diagram  in  Figure  18. 

For  calculation  of  the  total  missile  kill  probability  (  TPMK  ): 

PMK,  -  (I  -  r!*)(  (|  fDQD}  )PHK  , 

PMK,  =  (1  -  <T2“)(  (j  ^ D  )(1  -  PHK)PSPL  , 

PMK,  -  (1  -  e-2l)(  )d  -  (1  -  B)m)PHK  , 

PMK,  =  (l  -  <T2’)(  )(1  -  (1  -  «)M5)(1  -  PHK)PSPL  , 

PMK,  =  (1  -  *-*“>(-,  +  opf  )d  -  B)M5PHK  , 

PMK*  »  e~2*PHK  , 

TPMK  =  PMKX  -f  PMK2  +  PMK3  +  PMKA  +  PMK5  +  PMK* 

TPMK-(  I  -  ,r2*)(  (1  QDqd)  )PHK  +  (1  -  r2*)(  j-  )(l  -  PHK)PSPL 

+  (l  -  e'”><  TTTW  »(1  -  "  -  + "  -  TiTW  > 
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Figure  18.  Probability  calculations  for  one  line  ol  encounters. 


(1  -  (1  -  *)W5)(1  -  PHK)PSPL  +  (1  -  e~2*)(  ■  +'QD)  )(1  -  B)msPHK 


+  ( e~2*)PHK  . 

For  calculation  of  the  total  missile  diversion  probability  (  TPDIV  ): 

PDIVX  =  (1  -  e~2a)(  '{{~DqD)  )(1  ~  PHK)([  -  PSPL)  , 

PDIV2  -  (1  -  e'2*)(  —  )(1  -  (1  -  S)M5)(1  -  PSPL)  , 

PDIV3  =  e~2\l  -  PHK)  , 

TPDIV  -  PDIVX  +  PDIV2  +  PDIVi  . 

TPDIV- (1  -  e~2*){  ^Dq5)  )0  -  PHK){\  -  PSPL)  +  (1  -  e~2a)( )  • 

(1  -  (1  -  5)ws)(l  -  PSPL)  +  (e"2a)(l  -  PHK)  . 

For  calculation  of  the  ship  hit  probability  and  ship  survival  probability: 

PSHI  =  (1  -  <r2a)(  1  )(1  -  5)ws(l  -  PHK)  , 

PSSUR  =  (1  —  (1  —  e“2s)(  Jf+Q-D)  )(1  ”  B)MS{1  ~  PHK))  • 
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